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We investigate superconducting pairing driven by electron-electron interactions in a theoretical 
model for monolayer M0S2 with the temperature-flow functional renormalization group(fRG). At 
low doping, the dominant instability is toward odd-parity pairing with /-wave Mo-nearest-neighbor 
structure. We compute the fRG phase diagram versus electron doping below the van Hove filling of 
the conduction band. In the superconducting regime, the critical temperature grows with doping, 
comparable to the experiments. Near van Hove filling the system favors a ferromagnetic state. We 
demonstrate that the triplet pairing is driven by ferromagnetic fluctuations and that the multiorbital 
nature of the conduction band as well as the Hund’s coupling appear crucial in making the physics 
of M0S2 different from e.g. doped graphene. 
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Superconductivity with critical temperatures around 
lOK has been reported in M0S2 flakes by electron doping 
using combined ionic-liquid and solid-state gating (ij, [3. 
The observed superconducting dome as function of dop¬ 
ing suggests the possibility of unconventional pairing. As 
Mo is a transition metal and the electronic bands near the 
Fermi level in M0S2 have predominant d-orbital weight, 
electron-electron interactions are likely to play a role. 

Monolayer M0S2 in the 2 H-structure can be regarded 
as a sibling of of graphene. Compared to graphene with 
its two C sublattices, in M0S2, one sublattice site is occu¬ 
pied by Mo, while the other sublattice is formed by two S 
atoms displaced above and below the Mo plane. The lay¬ 
ered 2 H-structure is a direct gap semiconductor with gap 
size ^ 1 . 8 eV[l, 0 . Bulk M0S2 has an indirect gapjl^. 
Upon electron doping, the Mo-d-dominated conduction 
band develops pockets around the K and iA'-points. By 
continued doping, these pockets merge with small extra 
pockets around the M-points. At the merging points, the 
density of states (DOS) exhibits van-Hove singularities. 

The possible nature of superconductivity in M0S2 is 
of high interest to theory [j-Q. The work by Roldan et 
al.Q compares pairing due to electron-electron interac¬ 
tions and mediated by phonons using effective intra- and 
inter-pocket interactions between the pockets that form 
in the conductions band around the K and K' points 
upon doping. Their conclusion is that phonon-mediated 
interactions have a smaller effect, and that the resulting 
pairing is driven by the short-ranged Coulomb interac¬ 
tions, and comes out as odd parity. Topological p-wave 
states due to electronic interactions are also reported by 
a mean-field theory Q. A purely phonon-mediated mech¬ 
anism based on first principles is proposed in Refs. | 3 ^- 

The purpose of this paper is extend the study of Ref. 
H and to explore possible ground states of M0S2 using 
a purely electronic model. We apply a temperature-flow 
functional renormalization group (fRG) technique [l^. 
to a three-orbital tight-binding model for M0S2 proposed 
in Ref. [^. The T-flow fRG is an unbiased approach 



FIG. 1 : Left: Band structure along the contour F —A —M — F. 
Inset: Schematic view of crystal structure of monolayer M0S2. 
The purple sphere represents the Mo, the yellow spheres sul¬ 
fur atoms. Right: Electronic density of states (DOS) of the 
3 -band model. The red line shows the total DOS with the 
van-Hove singularity at A = 2 . 2 eV. The blue line is the 
band-projected DOS for the dUorbital, while d„,2_y2- and 
da;;/-contributions are shown by the coinciding solid purple 
and green dashed lines, respectively. 


and sufficiently flexible to treat the doping-dependent 
Fermi surfaces with enough angular resolution such to 
disentangle different pairing types. Importantly in our 
context, it is not blind with respect to ferromagnetic 
fluctuations [l^ . Furthermore, we use as input effective 
interaction parameters obtained by constrained random- 
phase-approximation (cRPA) calculations. Hence, in ad¬ 
dition to confirming the qualitative picture for the paired 
state of Ref. 0 , we give an ab-initio-based estimate of the 
critical temperatures as function of the doping, unveil the 
origin of electronic pairing as mediated by ferromagnetic 
fluctuations and the importance of Hund’s rule coupling, 
and give an overview of the competing phases. This 
suggests experimental searches for ferromagnetic corre¬ 
lations in M0S2 and related systems. 

For the low-energy physics in Mo^ we use a three- 
band tight-binding model from Ref. l 9 |. In the localized 
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state basis ^ , d^j. 2 _y 2 , ^ the tight-binding model 

Hamiltonian is given by a 3 x 3 matrix. For simplicity, 
we only consider nearest-neighboring hopping terms with 
the parameters of Ref. @ , resulting in bands as plotted 
in Fig. 1(a). Spin-orbit coupling will be discussed more 
below. Doped electrons enter first into the middle band 
around K(K'), forming Fermi pockets. In Fig.1(b), we 
plot the density of states (DOS) of this band. It exhibits 
a van Hove (VH) singularity around E = 2.21eV, which 
is when the extending Fermi pockets touch developing 
small pockets around M. In this paper, we focus on 
smaller doping levels x up to van-Hove filling. Around 
K, K\ the doped electrons have predominant charac¬ 
ter, while near the VH level, the and dxy compo¬ 

nents are stronger, as can be seen from Fig. 1(b). In the 
DFT band structures and more elaborate tight-binding 
parametrizations[^, additional band minima along the 
T-K direction not captured in our simplified model can 
lead to additional Fermi pockets at higher x. These how¬ 
ever have smaller DOS than e.g. the M-point regions and 
should therefore be of minor importance for the trends 
found in this work. The interactions are given by 

EE “1“ ^00 EE 

i a i a<b 

+ Jh E E ^ias fifes' 

i a<.b,s,s' 

+ E E + H.C^ 

i a<b 


with Uqo = 2.61eV, Uqq = 2.09eV, Jn = 0.27eV, com¬ 
puted from first principles by cRPA 12h14| . In a broader 
view, besides the minor splitting of the van Hove singu¬ 
larities near the M-points, the Fermi surface evolution 
in our model for M 0 S 2 is reminiscent of that in doped 
graphene (e.g. studied in [l9l - Elj |). The natural question 
arises if the phenomenology is similar. We we argue that 
the multi-orbital nature due to the mixing of 3 d-orbitals 
at the Fermi level in the conduction band causes qualita¬ 
tive differences, at least for model parameters close to the 
ab-initio predictions. We in fact show that the Hund’s 
coupling is a tuning parameter that takes us from the 
’doped-graphene’ scenario to M 0 S 2 . 

We use the temperature-flow fRG [l5l . [l^ to analyze the 
possible instabilities due to interactions. The main idea 
of the method is that at very high temperatures T the in¬ 
teractions of the model are given by the bare ones, as all 
loop corrections are killed by T. Upon lowering T, the T- 
flow tracks the build-up of the loop corrections. The main 
object of study is an interaction function VT(ki,k 2 ,k^) 
that depends in our case on three wavevectors ki , ^2 and 
k^ for two incoming legs 1 and 2 and one outgoing leg 
3 (with the same spin projection as leg 1) in the con¬ 
duction band. The flow equations for Vt(^i, ^ 2 , ^ 3 ) for 
charge and spin-rotational invariance are taken from the 


general fRG formalism for one-particle irreducible vertex 
functions (for a review, see Refs. EH or also Ref. 
[3) and are given in the supplementary material. We 
start the flow at an initial T of roughly twenty times the 
bandwidth. Lowering T generally leads to a strong dif¬ 
ferentiation of the components of VT(fci 7 ^ 2 , ^ 3 ), where 
certain components become large at low T and actually 
diverge. These generalized Gooper instabilities can then 
be interpreted as ordering tendency in a particular chan¬ 
nel, depending on which combination of fci, ^2 etc. ac¬ 
tually grows most strongly. The flow is stopped when 
the maximal component of Vr grows larger than thirty 
times larger than the bandwidth. This dehnes a critical 
temperature Tc which should be understood as upper es¬ 
timate for any true ordering. More information about the 
interpretation of the fRG flow can e.g. be found in the 
reviews of Refs, [l^ and [I 3 . In a true two-dimensional 
system, collective effects should in fact destroy any long 
range order that breaks continuous symmetries. These 
effects are however not included in the used approxima¬ 
tion in the fRG, so do not play a role in our results. Also 
in view of the embedding of the experimental systems 
in three dimensions we do not expect that such collec¬ 
tive effects have to be considered at this stage. Further, 
we patch Brillouin zone around the K, K' points with 
patch center points numbered by r = I..., V, as seen 
in Fig. [21 Then we discretize the coupling function as 
Rr(D,T' 2 ,T- 3 ), i.e. with components. At each corner 
K, K', we take 5 patches (schemes with more patches 
have also been checked, the result is stable). We divide 
the radial direction of each patch into three sectors, pro¬ 
viding some more resolution in the radial direction. In 
the loop calculations, we use 3 radial integration lines in 
each patch, but did also checks with more lines. 

We run the T-flow fRG for a wider range of elec¬ 
tron dopings X. For small x, we find pairing instabili¬ 
ties with total-momentum-zero Gooper pairs in the odd- 
parity channel, while at larger x near van Hove filling, a 
ferromagnetic (FM) instability is found. The latter can 
be understood as a renormalized version of a Stoner in¬ 
stability due to a large DOS at the Fermi level. 

The data for a prototypical pairing instability at x = 
0.03 is shown in Fig. |2]on the right side. The inset shows 
the coupling function for incoming and outgoing legs at 
the Fermi surface, as function of the two incoming patch 
indices, with the first outgoing leg fixed in patch 2. The 
two diagonal features are exactly on the lines where the 
two incoming wavevectors add up to zero. These pair¬ 
scattering processes grow to very large absolute values 
and cause the instability. We also see the sign struc¬ 
ture with pieces of large positive and large negative val¬ 
ues. This is detailed more in the main part of Fig. |2l 
right side, which is basically a cut along the zero-total- 
momentum diagonal features in the inset. By looking up 
the patch numbers Fig. 1 one extracts that the pairing 
changes sign npon rotating the incoming pair by 60 de- 
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FIG. 2: Left: Patching scheme used in the £RG. At low doping 
X the Fermi surface (FS) is around the K, K' points. On the 
FS we define 30 patch points (grey filled circles). The patch 
index r runs from 1 to 30. In addition to these FS points, we 
define two other sets of patch points (blue and pinks circles). 
This way, the energy dispersion is monotonic in each sector 
in radial direction. Right: Pair scattering at small x. The 
green symbols in the main plot show the fRG data for the 
low-T pair scattering V{t, f, 2) near the instability, where the 
two incoming wavevectors belonging to index r and f add up 
to total zero momentum. The first outgoing wavevector is 
fixed at patch point 2. The red square symbols show a fit 
to the fRG data using a nearest-neighbor-(i.e. Mo-Mo-)bond 
pairing —/p(k[r])/p (k[2]). Inset: P(ri, r 2 , ts) as function of 
patch indices ri and T 2 for rs = 2. Clearly visible is that 
only interactions with a special relation between ri and T 2 
(implying zero total incoming wavevector) become large. For 
the data in this plot, fj, — l.SeV. 


grees, from one Fermi pocket around, e.g. around K, to 
another pocket, e.g. around K'. The same sign change 
is found when the incoming pair is fixed and the out¬ 
going pair is varied. The couplings do not vary much 
around K(K') points, i.e. is no nodal structure along 
the Fermi pockets. This odd-parity pairing can be iden¬ 
tified to transform according to the one-dimensional A 2 
irreducible representation of point group of M 0 S 2 . In 
real space it can be understood as having its main con¬ 
tributions on the bonds from one Mo site to its 6 nearest 
neighbor Mo sites, with the sign alternating upon 60 de¬ 
grees rotation, implying the odd parity. The form factor 
in wave vector space can also be interpreted as /-wave 
and reads fp(k) = sin [VSkx) — 2 sin ^ ^ cos 

The ht to this form factor shows good agreement (red 
squares and green circles in Fig.3). 


We also run a momentum-cutoff fRG flow[l^ which 
disfavors FM fluctuations with small wave-vector trans¬ 
fer across the Fermi surface down to very low scales 151. 
These flows do not give any instabilities but show some 
tendency to even-parity pairing. Therefore, we can con¬ 
clude that /-wave pairing is driven by Q = 0 ferromag¬ 
netic fluctuations. A similar conclusion is reached by 
turning off specihc channel spin fluctuation channels in 
the flow equation directly. 

In Fig. Owe show the fRG phase diagram of our M 0 S 2 
model with doping. In the low doping regime, the dis- 


FIG. 3: Left: £RG phase diagram of monolayer M 0 S 2 versus 
electron doping x. The blue dots give the TcS according to 
fRG for the /-wave superconducting phase (SC). The pur¬ 
ple dots represent the TcS ferromagnetic state (FM), now 
scaled down by a factor of 100. The fRG T^s should be in¬ 
terpreted as upper bounds for true ordering temperatures. 
Right: fRG critical temperatures as function of Hund’s cou¬ 
pling Jh for both the triplet pairing (blue circles,) and FM 
instability (pink circles, again multiplied by 0.01). Notably, 
for Jh ~ 0 the dominant pairing instability appears in the 
d-wave singlet pairing channel. 


cussed /-wave superconductivity is the dominant insta¬ 
bility, with a relatively low critical temperature Tc (of 
order l.OmeV). We can also see that the critical temper¬ 
ature grows with increasing electron doping. The highest 
scale for pairing for our parameters is around S.OmeV, 
considerably higher than the experimental data around 
1.3meV. We comment on this further below. The op¬ 
timal doping for the /-wave pairing is reached around 
X ^ 0.1. This shows good agreement with experiment [ij. 
The increasing DOS with growing x also enhances the 
FM fluctuations. In simple mean-field or Stoner pictures 
the critical interaction strength required for FM order 
is inversely proportional to the DOS at the Fermi level, 
and as the DOS diverges at the van Hove filling, one 
will definitely find ordering in these approaches. In the 
fRG calculations we also obtain a regime where the lead¬ 
ing instability is of FM type, however in reduced extent 
compared to mean-field theory due to the competition 
with other fluctuations. Hence, the prediction form the 
fRG approach would be that there are predominant FM 
fluctuations around the van Hove filling. Whether the 
ordering actually occurs or not cannot be answered by 
our calculations. For the Hubbard model on the square 
lattice at large second-nearest-neighbor hopping, the fRG 
also predicts FM ordering, in agreement with variational 
Montecarlo results [ 2 ^. However, the van Hove filling 
may be hard to realize precisely in the experiment as dis¬ 
order and spatial fluctuations will tend to smear out the 
diverging DOS. Nevertheless, the FM correlations should 
still be visibly enhanced. 

At small X the states near the Fermi level are mainly 
of dz^-type. Thus the local intraorbital repulsion, Uqo, 
should be important. We hence run the flow with only 
Uqo being nonzero. As expected, the flow still ends up 
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with the /-wave instability. However, is lowered from 
about 4.8meV with full interactions down to 1.2meV. 
This shows that interorbital processes contribute signif¬ 
icantly to the pairing. Next we vary the Hund’s cou¬ 
pling Jh in the interactions as given in Eq. (3), with 
C/qo and C/qq fixed, and the doping fixed at /i = 1.8eV. 
From Fig. 4, we can read off an enhancement of Tc 
in the /-wave regime by a larger Jh, from S.OmeV to 
around S.OmeV. Thus, Jr enhances the /-wave insta¬ 
bility as well. Even more strikingly, when we set the 
Jh = O.OeV, the system gives way to an even-parity 
pairing structure, with most attractive eigenvalues in the 
d-wave symmetry channel [l^. associated with the two- 
dimensional E' irreducible representation of D^h point 
group. Complex combinations of dxy and break 

time-reversal symmetry, leading to chiral ed ge m odes as 
known from strongly doped graphene models [l9|,l23j. We 
also study the importance of the Hund’s coupling Jr in 
the FM regime. The data is shown as pink dots in Fig. 5. 
The critical temperature is again enhanced dramatically 
when Jh is increased. The Hund’s coupling intensifies 
the FM fluctuations just as it strengthens the /-wave 
pairing. Such an effect of Jr also reported in models for 
NaCoOa 2/^2 0 [ill. We also observe a kink in the critical 
temperature as a function of Jr at Jr ~ O.OSeV where 
the system undergoes a phase transition from even-parity 
pairing at small Jr to odd-parity pairing. Hence Jr has 
a key role as a switch between two different situations. 
On the one hand, we have the single-orbital scenario of 
graphene, where Jr is ineffective as the bands near the 
Fermi level are just made from C pj-states. Then one 
finds strong antiferromagnetic fluctuations and chiral d- 
wave singlet pairing. For M 0 S 2 , the multi-orbital char¬ 
acter of the conduction band makes Jr important and, 
for realistic values of Jr according to cRPA, one finds 
predominant FM fluctuations and triplet pairing. 

These considerations indicate that the superconduct¬ 
ing ground state is quite susceptible to parameter 
changes. This is one reason why the rather large pairing 
scales ~ lOmeV exceeding the experimental scales by one 
order of magnitude should not be taken as a failure of this 
theoretical picture. Some reduction of the TcS should be 
expected if we would be able to include self-energy correc¬ 
tions or more collective fluctuations into the fRG, which 
is however not yet technically feasible for multiband sys¬ 
tems. In view of the importance of the Hund’s rule cou¬ 
pling, we refer to recent works by Georges et al. , who 
pointed out that Jr typically degrades the quasiparticle 
weight resulting in a bad Hund’s metal. This may be an¬ 
other important factor that reduces the superconducting 
TcS compared to our study. Then there are uncertainties 
about the actual experimental system. Disorder is ex¬ 
pected to have an effect on the Tc of any unconventional 
superconductor, although in the /-wave case only large- 
momentum impurity scattering will have an impact of 
the phase structure of the pairing. On the positive side. 


the high scales from the theory might actually show that 
there may be some headroom for reaching higher T^s in 
Mo 5'2 and related systems. 

In M 0 S 2 and WS 2 , the^in-orbit coupling splitting can 


tne^u 

reach the order O.leV [1,0. In the conduction band near 
K, K', the actual values may be smaller ^ 3meV[^. 
The A^-component is conserved around the small Fermi 
pockets shown in Fig. 2 which then split into two pock¬ 
ets around each K, K'. By time-reversal symmetry, the 
pockets with the same radius at different pockets will 
have opposite spin projections. The f-wave pairing does 
not violate time-reversal symmetry 2^ and should not be 
changed strongly by the small splitting near K, K'. In 
such a case, the 2 ; component of the gap function d^(k) 
can be expected to be dominant. This corresponds to 
Sz = 0 pairing. 

In conclusion, we investigate the possible pairing struc¬ 
ture in monolayer-MoS 2 by the unbiased T-flow fRG. For 
interaction parameters obtained by cRPA, the leading in¬ 
stability at low electron doping lies in the /-wave pairing 
channel. We argued that ferromagnetic fluctuations pro¬ 
vide the pairing glue. The phase diagram exhibits dom¬ 
inant /-wave pairing up to x ~ O.IO. For doping close 
to the van-Hove filling, sufficiently strong interactions re¬ 
sult in FM order. The pairing Tc is maximal around an 
optimal doping of a; ^ 0.1, in quite a good agreement 
with experiment. Varying the Hund coupling strength 
Jr showed strong effects on both the /-wave pairing scale 
and the FM tendencies. Notably, when Jr is relatively 
small or set to zero, the leading pairing instability is in 
the singlet d-wave channel. This emphasizes the impor¬ 
tance of the Hund’s coupling and the multi-orbital nature 
of the Mo-d-dominated conduction band. 

Our study extends earlier work by Roldan et al.0 who 
already found /-wave pairing. It clarifies the nature of 
the pairing mechanism, the influence of the model param¬ 
eters, the importance of FM fluctuations and the Hund’s 
coupling. Future studies should include the phonon- 
mediated electronic interactions. This should allow one 
to compare the relative importances of these two distinct 
pictures for the superconductivity in M 0 S 2 . 
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Supplementary Material 
T-flow equations 


^ph,T{ki,k2,k3,k4) = - j dpL{p,p + p2 -pa) 
VTiki,p + k2 - k3,p)VTip, k2,k3) 


( 5 ) 


The temperature-flow fRG uses the physical temper¬ 
ature of the system as flow parameter. As in previ¬ 
ous works using this method the flow of the self¬ 

energy is neglected and the flow equations are truncated 
after the four-point vertex. Furthermore, the interaction 
vertex is approximated to be frequency-independent, and 
the external frequencies are set to zero. The wavevector 
dependence is treated in the A^-patch discretizationfisl- 

The initial condition for the flow at the initial T larger 
than the bandwidth of the system is the bare action, as all 
perturbative corrections vanish at high enough T. Then 
the flow leads to lower T and picks up the different one- 
loop corrections for the interactions, each in second order 
in the vertices at the given temperature. More precisely, 
the flow equations read in this casefl^ 


The ki indices are patch numbers, corresponding to 
wavevectors in the conduction band. In the tempera¬ 
ture flow without self-energy corrections, the one-loop 
integrals are given by 


L{p,p') = 


dT 


ILJ — Cr, liO — €r 


( 6 ) 




Y\d 

^ ph,T 


- ph,T 


where explicitly 


^pp,T{ki,k2,k3,k4) = - j dpL{p,ki + k2 - p) 
X VT{ki,k2,p)VT(j>, -p + Pi +P2,P3) 


•pd 
^ ph 


( 2 ) 


( 3 ) 


^T{ki,k2,k3,k4) = - j dpL{p,p + pi -pa) 

X [-2VTiki,p, k3)VT{p + ki- k3, k2,p) (4) 

-I- VT{ki,p,p + ki - k3)VT{p + ki - k3,k2,p) 

+ VT{ki,p, fc3)VT(fc2,P + fci - k3,p)] 


In most situations, the usual momentum-shell cutoff 
schemes which have a cutoff-derivative in the loop di¬ 
agrams give equivalent flows as the TfRG. In the case of 
particle-hole fluctuations with small wavevector transfer, 
a momentum-shell cutoff would lead to underestimation 
of these processes at intermediate scales while the T-flow 
incorporates these contributions in the same way as fluc¬ 
tuations with larger wavevector transfers. Hence the T- 
flow can capture ferromagnetic instabilities and is hence 
a more unbiased fRG-method. Note that also frequency 
cutoff can be used to circumvent this shortcoming of the 
momentum-shell schemes 21. MlIS- 









